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Abstract 

We consider packings of congruent circles on a square flat torus, i.e., periodic (w.r.t. 
a square lattice) planar circle packings, with the maximal circle radius. This problem 
is interesting due to a practical reason the problem of ''super resolution of images." 
Wc have found optimal arrangements for = 6, 7 and 8 circles. Surprisingly, for the 
case N = 7 there are three different optimal arrangements. Our proof is based on a 
computer enumeration of toroidal irreducible contact graphs. 

1 Preface 
1.1 History 

A class of problems of an optimal (in different senses) placing of equal balls in different 
containers is a huge class of interesting research problems in discrete geometry. The best 
known variations are the following: to find the densest packing of equal balls in M", to find 
the minimal covering of by equal balls, to find the maximal number of equal balls that 
can touch another one with the same radius, and to find the maximal nunAer of equal 
balls that can fit into a prescribed container. One can also allow balls of different radii or 
even replace the balls with something hke simphces or other objects, or study the universal 
properties of packings of convex bodies, giving rise to many other similar problems that are 
also of interest. These problems have been widely studied and appeared to be very complex 
despite having very simple formulations. Another point is that many of these problems have 

*The author is supported by the Delone Laboratory of Discrete and Computational Geometry (P. G. Demi- 
dov Yaroslavl State University) under RF Government grant 11.G34.31.0053, RFBR grant 11-01-00735 and 
NSF grant DMS-1101688. 

^The author is supported by the Chebyshev Laboratory (Department of Mathematics and Mechanics, St. 
Petersburg State University) under RF Government grant 11.G34.31.0026 and by the Delone Laboratory of 
Discrete and Computational Geometry (P. G. Demidov Yaroslavl State University) under RF Government 
grant 11.G34.31.0053 



practical value (even for the "strange" containers and higher dimensions). Some of the results 
are collected in the following books: [T], [6], [10]. 

In this paper we are talking about one of the aforementioned problems, where the con- 
tainer is the flat torus /I? with the induced metric and the goal is to find the maximal 
radius of the disk, whose congruent copies can be fit in the torus without overlapping. 
The problem is equivalent to finding the densest packing, as the area of the torus is constant 
and equals one, and so the best density is 

P = NnRl^^^. 

We found out about this problem from Daniel Usikov ([19j), who described the practical 
reason for studying optimal disk packings on the torus. According to him, this case is espe- 
cially interesting for the problem of "super resolution of images." Here are several extracts 
from his letter that provide a more detailed description: 

Let's assume that the image obtained by an observation has an infinite ac- 
curacy. Now assume that we have just two observations. If both of them are 
conducted from the same point, obviously, these two observations are not better 
than just one. 

Now let's shift the position of the second with respect to the first one. Then any 
shift would be equivalent from the point of view of super resolution restoration. But 
in practical applications we need to take into account two sources of observational 
errors. First, the exact amount of the shift can be found with some accuracy only, 
and second, the brightness in each pixel of the image has an intrinsic error - so 
called "photon noise". So if the shift is within its statistical error, any attempt 
to calculate the super-resolution image will fail. Furthermore, it is visible from 
the formulas that they operate with the differences in brightness between adjacent 
pixels. So, by separating the observation points as far as possible, one can reduce 
the error in the estimation of the brightness differences. 

It definitely looks like the best result of super resolution observation would be 
to maximize the distance between two points of observations. For the case of the 
square pixel, it is the half of the pixel's diagonal. When a third point is added, we 
expect that the best allocation would be to maximize the minimal distance between 
three point, etc. 

I think that a rigorous prove of the max-min requirement is achievable. As 
a simplification on the road to the prove, one can assume that the brightness in 
each pixel takes a random value within an interval (say 0:1). Then averaging 
over all possible random brightness samples, the brightness should gone from the 
optimization formulas, and the task will be reduced to the search of the optimal 
location of observation points on the torus. 

Sent by email, from Daniel Usikov (dusikov@yahoo.com) 
November 22, 2010 
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Obviously, the optimal packing in the torus could not be worse than the optimal packing 
in the square [0, 1]^. Here are some results for the small number of disks in the square, 
borrowed from [1] (See [12], [T3], |15j). Here d denotes the distance between the centers. 
Corresponding configurations are shown in Figure [TJ 
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Figure 1: The optimal configurations for the square 



1.2 The formal statement 

The problem is: for a given number > 1 of points, find the maximal r G such that 
circles of radius r could be put on the square fiat torus T = M^/Z^ without overlapping, 
or, equivalently, to find the maximal d G such that there are points on the torus with 
pairwise distances not less than d (where d = 2r). 
We will work with the second statement. 

1.3 Known results for the torus 

Optimal configurations of up to 4 circles are relatively easy to find, and for < 5 they 
were found in [8j. These configurations and the corresponding d-s are presented in Figure [2] 
D. Usikov (p,9j) has conjectured some optimal configurations using numerical simulations 
for the small numbers of disks. 

2 Main results 

In this paper, we present solutions of the problem for A^ = 6, 7 and 8 and a conjecture for 
the case A^ = 9. 
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iV = 2 /V = 3 r , ^ = 4 (A, ^ = ^ 
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Figure 2: Optimal configurations of up to five points in tlie torus 



Theorem 1. The arrangement of 6 points in which is shown in Figure 3a is the best 
possible, the maximal arrangement is unique up to isometry, and d ^ 0.40040554. The precise 
formula is 

VS 72^373 + 2 1 
~2 
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Theorem 2. There are three, up to isometry or up to a move of a free disk, optimal ar- 



rangements of 7 points in which are shown in Figures 



3b 



3d 



where d 



0.3660. 



Theorem 3. There is one unique, up to isometry, optimal arrangement of 8 points in T^, 

1 



which is shown in Figure 



3e 



where d 



0.3660. 



Conjecture 2.1. There is one unique, up to isometry, optimal arrangement of 9 points in 
T^, which is shown in the Figure 



3j where d 



0.3437. 



\/5+2V3 

So, for = 6, 7 and 8 we have a 5 — 7%-improvement compared to the results for the 
square. The exact point coordinates can be found in the Appendix [B} 



3 Basic definitions and preliminary results 

The main notion we use to find an optimal configuration is the notion of a contact graph. 

Contact graphs. Let X be a finite set in T^. The contact graph CG(X) is an embedded 
graph with vertices X and geodesic edges (x, y), x,y ^ X such that x and y are adjacent iff 
dist(x, = min dist(M,f). 

As a corollary of the definition, one sees that all the edges in a contact graph have equal 
lengths. 

We refer to the configuration of A^ points that maximizes the minimal distance between 
the A^ points, i.e., solves the packing problem, as optimal. To avoid dealing with an infinite 
family of torus isometrics induced by parallel transports of the covering plane, from now on 
we suppose that every point set in question contains the vertex (0, 0). 
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(b) iV = 7, optimal 




(d) N = 7, optimal 



(e) = 8, optimal 



(f) = 9, conjecture 



Figure 3: Configurations for = 6, 7, 8 and = 9 
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Irreducible contact graphs. We say that a contact graph CG(X) (and its corresponding 
configuration X) is irreducible if it is rigid in a sense that there are no possible slight motions 
of vertices (we do not distinguish the points in X and the vertices of the graph) that change 
the contact graph (either the graph itself or the topology of its embedding) or increase the 
edge lengths. A more detailed definition of the contact graph for the torus using so called 
strut frameworks and examples can be found in [8] (there it is called "infinitesimaly rigid 
strut framework"), which as they state is equivalent in our case to the definition in [3j. The 
only difference between our definitions is that we allow an irreducible contact graph to have 
isolated vertices. 

This definition can be applied to point configurations in any metric space. Note that 
W. Habicht, K. Schiitte, B.L. van der Waerden, and L. Danzer used irreducible contact 
graphs for the kissing number and Tammes problems for [TT| [T6| IT7\ [7j. The idea of 
irreducibility is the expression of the thought that the optimal configuration should have 
some nice properties and there is a possibility to "adjust" any optimal configuration to have 
even more nice properties. Here we only want the configuration to be rigid, but, for example, 
Danzer [7] used a trick to disallow the optimal graph to admit another adjustment that is 
now called "Danzer fiip", and that helped him to solve the Tammes problem for 10 and 11 
points. Also this trick helped to kill lots of graphs in [13j . 

Now we state some basic propositions that will be of significant use. 

Proposition 3.1. Let X be an optimal (corr. locally optimal, i.e., giving the local maximum 
of d in the natural topology) configuration of N points on a square flat torus. Then its contact 
graph is either irreducible, or X can be continuously moved to another optimal ( corr. locally 
optimal) configuration with an irreducible contact graph. 

Although this fact is easy to prove, we refer to [3]. Also see Theorem 3.2 in [8J. 

Proposition 3.2. For any configuration X its contact graph can have vertices of degree at 
most six. 

Proof. It is obvious since any circle in the plane can not touch more than 6 circles of the 
same radius, and since the universal cover of the contact graph is the contact graph in the 
plane, the property holds. □ 

Proposition 3.3. For any configuration X its contact graph can not have angles between 
adjacent edges below |. 

Proof. It is obvious since an angle must be not less than in an equilateral triangle. □ 

Proposition 3.4. There are no angles greater or equal than it in an irreducible configuration. 

Proof. If there is an angle of tt or greater, than the corresponding vertex can be slightly 
"pushed" inside a face with that angle, and it becomes isolated. □ 

Corollary 3.1. There are no vertices of degree one or two in an irreducible configuration. 
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Proposition 3.5. An irreducible contact graph without isolated vertices is cellular (i.e., the 
interior of each face is homeomorphic to an open disk). 



Proof. If not, then one can easily see that such graph has an obtuse or a flat angle. For the 
details see fS|. □ 

Proposition 3.6. An isolated vertex in the contact graph can be only inside a face, whose 
universal cover is a polygon of at least seven vertices. (The universal cover is referred to 
for the strictness, as there can be, for example, "octagonal" faces on the torus with only 4 
distinct vertices. Nevertheless, in the following we will call such faces just "octagons".) 

Proof. If we draw all straight-line segments from the isolated vertex v to all the vertices of 
the face, all of the segments will be strictly longer than the edges, thus the angles are strictly 
less than |. Their sum is 27r, so there must be at least 7 such angles. By the way, if there is 
a vertex inside a hexagon, then it is a regular hexagon split into 6 equilateral triangles. □ 

Remark. If we remove an isolated vertex from the contact graph on vertices, then we 
will get a locally optimal configuration on — 1 vertices. 

Proposition 3.7. An irreducible contact graph without isolated vertices has at least 2N — 1 
edges. 

Proof. We refer to the result of Connely |5]. Some general position ideas that can be 
used to prove this will be given in the next sections. □ 

Now we state some propositions involving the d we are looking for. 

Proposition 3.8. Let d{N) denote a solution to the packing problem for N points. Then 
d{N) is decreasing. 

Proof. One can remove any vertex to achieve a configuration on — 1 points. □ 

2 

Proposition 3.9. The following inequality holds: d{N) < 



/NV12 

Proof. It follows from the knowledge of the maximal density of circle packing on a plane, 

which is (p], [10]). See |8]. □ 

V 12 

Proposition 3.10. There are no multiple edges in the contact graph for N > 5 points. 

Proof. There are no multiple edges in the universal cover since any two disks on the plane 
cannot touch in two points. So the multiple edges Ci = (uv) and 62 = (vu) may occur only if 
the cycle uvu = 616261 wraps around the torus (i.e., is not homologically trivial). But then 
the length of the edge must be at least 1, which is not possible as d{5) = ^ < ^ (|H]) and 
d{N) decreases. □ 
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4 The idea of the proof 



Now we have some basic properties of contact graphs, and so we are going to use a computer 
to find the one that will give the optimal configuration, among all the graphs with such 
properties; the idea is similar to the one used for solving the Tammes problem in [13]. There, 
to find an optimal configuration on the sphere, all the planar graphs with some necessary 
properties have been enumerated and tried to be embedded into a sphere with equal edge 
lengths, satisfying the properties of contact graphs. The enumeration has been done by a 
program plantri by G. Brinkmann and B. McKay([2i|). 

We've used the program surftri written by T. Sulanke ([IH]), who very kindly modified 
the original release, that enumerated only the triangulations, to be able to enumerate all 
toroidal cellular graphs, and we are very grateful to him for that. 

But unlike |13| . where all planar graphs have been enumerated and the possibility of 
irreducible embedding has been checked, we've used a slightly different approach. The idea 
grows from the fact that since we are unable to find an embedding of the graph exactly, in 
any case it would be very hard to distinguish the graph and its subgraphs programmatically. 
Here are the details. 

Suppose we have an embedding of the graph G without isolated vertices, like in any of 



the figures (except for Figure 3b) above. Then, considering each edge as a vector, we may 
write the following system of equations that characterizes the embedding (up to a parallel 
transport that is out of interest for us): 

' \v,\^ = d^ i = l..E 

i^j = f is a face of G 

Vi is an edge of f 

±Vi = const 1 (1) 

Vi is an edge of ci 

^ ±Vi = const 2. 

^ Vi is an edge of C2 

Here we assign a vector Vi to each edge of G, fixing any of the two possible directions of 
the edge. Then the first E equations tell that all edge lengths are equal to d. The next F 
equations say that each face is really a face, i.e., the sum of the vectors corresponding to 
its edges (with appropriate signs) equals to zero. The last two equations "fix the cut" of 
the torus: we choose any two cycles Ci and C2 that do not form a face (i.e., homologically 
independent) and prescribe for each of them the way it wraps around the torus. Since it is a 
cycle, then the sum of the vectors, corresponding to its edges, has integer coordinates, and 
based on the number of the edges in the cycle, we easily establish that there is only a finite 
number of possible values of consti and const2 if the cycle lengths are bounded. The details 
will follow in lemmas [l}]4j An example of the system for a specific graph is System ^ on 
page[T7} 

Remark. It is well-known that a graph embedded on a surface is unique up to orientation- 
preserving homeomorphism of the surface characterized by its rotation system — a cyclic 
"clockwise" order in each vertex of the out-going edges. And surftri outputs these rotation 
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systems. However, there are many homeomorphisms of the torus from SL2(Z) that do not 
conserve lengths, and so we have to somehow "fix the homeomorphism", and in the previous 
system it is done by selecting the two exclusive cycles and fixing them. The details on that 
will be below. 

Note that the equations on the faces are not linearly independent (except the case where 
there is only one face, which is trivially impossible with our restrictions on degrees of vertices 
and without multiple edges): each edge is counted twice, once in positive direction and once 
in the negative, so the sum of all these equations is zero. Removing one of the equations 
for the faces, and, rewriting the system in coordinates, we get a system of m equations in 
n = 2E + 1 variables (coordinates and d), where 

m = E + 2{F - 1) + 2 ■ 2 = E + 2F + 2. 

By the Euler formula for a torus, we have a relation N + F = E, so 

m = E + 2{E - N) + 2 = 3E - 2N + 2. 

If a solution exists, by the reasons of general position, we may hope to be able to find it by only 
considering only nofm equations. On the geometrical level, it may be the same as considering 
the subgraph instead of the original, hoping to get the removed edges automatically. Taking 
m = for a subgraph we find that E = 2N — 1. By Proposition |3.7 the contact graph has 



at least such number of edges, so there exists its subgraph having exactly 2N — 1 edge^ So, 
our idea is to consider only the graphs with 2N — 1 edges and to find their embeddings. 
Remark. Although the system of equations is polynomial and there is much known about 
them, algebraic surfaces given by these equations seem not to be in the general position in 
the whole space, and such a system usually has infinite number of (complex) solutions, so 
the known approaches for polynomial systems fail. 



5 Technical details 

Now we need some statements to guarantee that the code for solving the problem works. 
Notation: For an (oriented) cycle c we write l{c) := Yl^i^ where Cj are the vectors Vi 
from System ([T|, corresponding to the edges of the cycle, with appropriate signs (that are 
determined by the orientation). 

Remark. In the following, we refer to a pair of (oriented) cycles, a cut along which unwraps 
the torus to a fiat polygon, as "homologically independent cycles", and to a cycle, that forms 
a face, as "trivial cycle". 

Lemma 1. If d < d(5) = then (in an irreducible contact graph) there is no non- 
trivial cycle of three edges where one edge occurs twice. Also, for any non-trivial cycle c of 

^By the way, the same reasons of general position help us think that if there are less then n equations, 
i.e., less than 2N — 1 edges, then the space of solutions has positive dimension. Thus there is an infinitesimal 
motion possible and a graph with less than 2N — 1 edges could not be an irreducible contact graph. 
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three edges l{c) = (±1,0) or l{c) = (0,±1), and for any non-trivial cycle s of four edges 

l{s) e {(±1, 0), (0, ±1), (1, ±1), (-1, ±1)}. 

Proof. It easy follows from such bound on d that these cycles lifted to just couldn't reach 
any other point in Z^: the distance in between the endpoints of any chain of three edges 
is not more than 3d, which is strictly less than \/2 — the length of a diagonal of the square, 
so, if this chain connects the two different lifts of the same vertex, then they could only be 
different by the unit horizontal or vertical vector. For a cycle of length 4, the argument is 
the same: the distance is not more than 4d < 2, so the shift is only possible by the unit 
horizontal or vertical vector or by the diagonal of the square with the length ^/2. □ 

Lemma 2. If d < ^ and Ci and are homologically independent cycles (in an irreducible 
contact graph) having lengths 3 and 3, then there exists an isometry of the torus (possibly 
orientation-reversing) that makes l{ci) — (1,0) and l{c2) — (0, 1). 

Proof. These sums are originally ±(1, 0) and ±(0, 1) or vice versa (by the previous lemma), 
since the reflections against the lines y — 0, y — x and x — are isometries of the torus, we 
easily achieve what we want. □ 

The next two lemmas have the same proof: 

Lemma 3. If d < ^ and c\ and ci are homologically independent cycles (in an irreducible 
contact graph) having lengths 3 and 4, then there exists an isometry of the torus (possibly 
orientation-reversing) that makes /(ci) = (1,0) and l{c2) — (0,1), or l{ci) — (1,0) and 
1{C2) = (1, 1), or /(ci) = (1, 0) and /(cs) = (-1, 1). 

Lemma 4. If d < ^ and ci and C2 are homologically independent cycles (in an irreducible 
contact graph) having lengths 4 o-nd 4, then there exists an isometry of the torus (possibly 
orientation-reversing) that makes /(ci) = (1,0) and /(C2) = (0,1), or l{ci) = (1,0) and 
1{C2) = (1, 1), or /(ci) = (1, 0) and /(ca) = (-1, 1), or /(ci) = (1, -1) and /(cs) = (1, 1). 

The next technical lemma guarantees some nice properties to the subgraph. 

Lemma 5. If N = 6 or N = 7, then for each irreducible (locally) optimal contact graph 
without isolated vertices there exists its subgraph with 2N — 1 edges which is cellular and all 
its vertices are of degree at least 3. 

Proof. We prove that if there is a cellular graph with degrees of all vertices at least 3 and 

at least 2N edges, then we can remove an edge so that the property holds. 
Firstly, we prove the estimate on the number of vertices of degree 3. Since 

3^3 + 4^4 + Sris + 6^6 = 2E, 

where is the number of vertices of degree k, and 

E>2N ^2n3 + 2n4 + 2^5 + 2^6, 
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we have 77.3 < + 2n6. 

Now we prove that there are at most two edges, the removal of which violates cellu- 
larity. Suppose there is such an edge e. Then there we have one of the following pictures 
(topologically ) : 




(a) 2 vertical (b) 1 vertical (c) 1 vertical 

1 ship 1 ship 



It means that there are several "vertical" components that wrap around the torus, several 
"ships" — planar components, and they are connected by the edges like e that cannot be 
removed without a loss of cellularity. 

Note that each "vertical" component may have at least three vertices as it should wrap 
around the torus and d < ^, and therefore at least four, as there are no multiple edges and 
all degrees are at least 3. So in the case of at least two different "vertical" components there 
should be at least 8 vertices, but we have only 7. If there is a "ship", then it has to contain 
at least 4 (even more, if we look more accurately) vertices since otherwise the third vertex 
has degree at most two. Since the embedding is cellular, there exists a "vertical" component, 
and thus we again have at least 8 vertices, which contradicts the statement. So there are 
no "ships" and in this direction there is at most one such edge (in a sense that if we cut 
the torus across this edge and get a cylinder then no other edge can be removed to make 
the embedding in the cylinder non-cellular). Similarly, inside the only 'Vertical" component, 
there is also at most one edge in the other direction, after removing which the embedding 
becomes non-cellular. 

If = 6, then we have ng = and 723 < 3 as 723 < ns and n-^ + 7i5 < 6. So we have at 
most 3773 < 9 edges prohibited to remove by the condition on degrees and at most 2 by the 
one on cellularity. So, among at least 2A'" = 12 edges, only 9-1-2 = 11 are disallowed to be 
removed. 

If = 7 and 776 >= 2 then, by the remark on the structure of vertices of the degree 6 the 
edge between these 2 vertices is a diagonal of a rhombus and can safely be removed. Huq = 1 
then we have 773-1-775 < 6 (the sum is less than 7 and must be even) and 773 < n^-\-2nQ = 775-I-2. 
Thus 773 < 4. Our graph has at least 2N — 14 edges, so in the case 773 < 3 we have at most 
3773 -I- 2 < 11 prohibited edges and we are done. The last case is 773 = 4. 
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If 77-6 = 1 and n3 = 4, then = 2 and these 2 vertices of degree 5 share a common 
adjacent vertex of degree 3, and all these 3 vertices are adjacent to the one of the degree 
6. That gives us that there is a rhombus abed, where dega = 5, deg6 = 6, degc = 5 and 
deg d = 3, from which we can safely remove the diagonal edge bd. That finishes the proof. □ 

6 Proof of Theorem [ll 

6.1 First step: finding all possible subgraphs 

To prove the desired Theorem[T| we ask surf tri to enumerate all cellular toroidal graphs on 
6 vertices having 6 ■ 2 — 1 = 11 edges and all vertices of degree no less than 3 (these are the 
subgraphs described in lemma |5]). There are 252 such graphs. Then for each graph we find 
two shortest homologically independent cycles and write down all the possible systems ([T| 
from page |8] (the possibilities are described in Lemmas [l] - 14]) and try to isolate its solution 
(or solutions) on the set 

\viy\ < d, i = 1..E 

d e [do, di] 
dist{pi,pj) > d, 

where do and di are the known bounds on d (the lower one comes from the "guess" for 
the optimal configuration, the upper one is theoretical). The isolation is performed using the 
branch and bounds method, i.e., on each iteration we subdivide the allowed intervals for each 
variable and check the possibility of the feasibility on the set using minimal and maximal 
estimates. We've used two different methods for choosing the independent variables to be 
localized by branching. 

In the first approach, to minimize the number of independent variables we first solve the 
linear part of the system concerning only x coordinates, where the space of the solutions has 
the dimension E — F = N, and we put bounds only on the coordinates Vj^,j G J that 
determine the solution (also, we put bounds on d), while the ^/-coordinates are determined 
from the equations Viy = i^/rf^ — fj^. The feasibility is checked in the following steps: 

1. Calculate the bounds on the Vi^ for i ^ J from the linear system on x-s; 

2. Calculate the bounds on the \viy\ corresponding to the already bounded x-s from the 
\vi\ = d equations; 

3. Enumerate all possible 2^ signs for Viy] 

4. Minimax check of the feasibility of the linear system on Viy. 

In the second approach, we solve the linear part on both x-s and y-s (i.e., on the vectors), 
and then use the expressions 

= d - cos(ai), 
Viy = d - sin(ai) 
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for the base vectors. And we put bounds on these Oj. Here the process of checking the 
feasibihty looks hke this: 



1. Calculate the bounds on the Vi^ and Viy for i ^ J from the linear system; 

2. Check lengths of these vectors. 

Remark. Actually, the second approach appears to work faster in general, since there is no 
need to enumerate all the possible combinations of signs for y-s, but, firstly, the problem for 
= 6 has been solved by using only the first approach, and, secondly, for a few graphs on 
7 vertices it also appeared to be faster. And the preliminary elimination for N = 7 has also 
been done using the first approach. 

If the solution is possible and the desired precision is not achieved, we continue the 
subdivision. 

In such a way, we find (subgraphs of) locally optimal configurations without isolated 
vertices with an arbitrary precision. Then, we do the same procedure for the graphs having 
less vertices and try to find some place for an isolated vertex there. 

This gives us some set of embeddings, but we are still unable to distinguish the subgraphs 
of the optimal configuration. So, at this stage the computation shows that the following 
theorem holds: 

Theorem 4. Each irreducible optimal configuration on the N = 6 vertices has a subgraph 
isometric to one of the graphs in the Figures^ and^ where each red square (that hardly 
could be seen on the image) determines the possible region for the point, and green circles 
are of radius about d. The precise numerical bounds could be found in the Appendix\A\ 




(d) Gi (e) G2 

Figure 4: Graphs pretending to be the subgraphs of the optimal 
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(a) Gs 



(b) G4 



Figure 5: Graphs pretending to be the subgraphs of the optimal 
6.2 Second step: distinguishing subgraphs 

To show that the only reason these graphs are found by the code is because they are the 



subgraphs of the optimal one showed on Figure 3a and do not lead to another optimal 
configuration, we distinguish two possibilities: if there surely is an obtuse angle in the graph, 
then the corresponding vertex should either be isolated (it is out of interest for us as in 
this case we will get this graph from the enumeration of graphs with less vertices), or there 
should be another edge from it. Looking at Figure [5} we easily see that for G3 and G4 there 
are obtuse angles, allowing the only possibility of adding an edge. If the edge is added, we 
obtain the graph shown in Figure |3a| where we can delete another edge from, and reduce 
the problem to, for example, the graph Gi. Now we are left with the two graphs Gi and 
To deal with them (as the pentagons are hard objects for manual maximization) we use the 
following techniques: we will prove that within the prescribed rectangles there is at most one 
solution of System ([T|, and since the one in Figure 3a satisfies these bounds, the solutions 



coincide. For convinience, we replace (P by just 2d, changing the bound in the appropriate 
way. 

To prove the propsition, first of all we rewrite System ([T]) in coordinates of the points of 
the configuration (we've multiplied by ^ for convinience): 



Xo = 

yo = (2) 
l{xi - Xj + + \{yi- Vj + lijf -d = 0, < i < j < N, {ij) is an edge 

where G {0,±1} shows which distance (of the nine possible) between the points we 

have to consider. We assume that the precision is good enough to determine these uniquely. 
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Suppose we have some 



-^0 — (3^1, Vi, 1/2, ■ ■ ■ , Xj^_i,yj^_i, d ) 

such that F{Xq) = 0, where F is the function M^^"^ — )■ M^^"^ on the left hand side of 
System g (we've substituted Xq = yo = 0) . Since F is quadratic and F{Xq) = 0, we have 

F{Xo + ev) = edFx,{v) + e^^d^Fx,{v, v) 

for any vector 

V = {dxi, dyi, dx2, dy2, ■ ■ ■ , dxN-i, dy^-i, dd) G B?^~^ . 
We claim that for e > small enough]^ 

||F(Xo + £t;)||oo>0 

for any unit (in || ■ ||oo norm) vector f , and if we would have an explicit bound on e then we 
could establish the uniqueness. 
Firstljj^ we compute 



and thus 



{d^F{v,v)){ij) = {dxi - dxjf + {dyi - dyjf 
1 



-d^F(v,v) 



< 4. 



Now, we have to estimate ||dF(f)||oo from below. 
We have 

(dFxo(f = (x- -x^j+ Sij){dxi - dxj) + {y^^ - + '^ij){dyi - dyj) - dd, 

where the coefficients are actually the coordinates of the edge vectors. We are interested in 

1 



inf \\dF{v) 
lklloo=l 



IdF^^II 

Xo Hoc 



(3) 



However, we do not know the coefficients exactly, so we need the following well-known 
lemma from the theory of linear operators: 

Lemma 6. For any two matrices A and B with A — B = A, the following inequality holds 

\B-m\A\\ 



\A 



-1 D-l 



5-^11 < 



1 - IIAIIII5- 



if the denominator is positive. 

^Recall that ||m||oo is the maximal of the absolute values of the coordinates of u. 
^F(^ij) means the coordinate that corresponds to the equation for the edge ij. 
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Applying this for A = dFxo and B = dFj^ where X is some point configuration, whose 
edges lie within the 5-neighborhood of that of Xq, we (noticing that since matrices A and B 
are different by 6 only in 4 positions in each row) get ||A||oo = ||y4 — 5||oo < 45 and hence 



|d-FY,!lloo < II -B """1100 + 



4:S\\B 



-1 1|2 
I oo 



xoiloo _ „^ ■ l-4(5||5-i||oo 

Substituting this into ([s]) we get the lower bound on dFxo(f); denote it by h. So, we have 

\\F{Xo + ev)\\ >eh-4e\ 

which is positive provided £ < |. So, given the bounding intervals for the edges we now can 
guarantee the uniqueness of the solution by calculating h and checking the inequalities. We 
do this for the graphs pretending to be optimal and Theorem [T] follows. 



7 Proof of theorems 2 and 3 



First of all, we notice that the Theorem|3] follows from the Theorem [2] because we know that 
d{8) < d{7) and there is a configuration with d = d{7). Uniqueness follows from the remark 
that after removing any vertex we would achieve an optimal configuration on 7 vertices, all 
of which are described by the Theorem [2| 

Now let's turn to = 7. We've used exactly the same method for finding the subgraphs 
as in the previous section. And (after a couple of weeks of calculations and code optimization) 



we've been left with the subgraphs of the ones in he figures 3b -3d Surprisingly, we had not 



been aware of the graph in the figure 3d, it has been found by the code. There are a bit 
more of the found subgraphs than for A^ = 6, so we do not provide all the images right here; 
nevertheless, they can be found at the URL below. There they are categorized by the way 
used to complete them to optimality and by the step of the elimination process. 

Now we have the last problem to distinguish the subgraphs. Note that the graph on 



Figure 3b is achieved by enumerating graphs on 6 vertices, and it has 12 edges, while we 
enumerate its subgraphs with 2-6 — 1 = 11 edges, so for it we need to prove the existence of 
the last edge. Similarly, we have to prove that the one last edge should be actually present 



for the subgraphs of [3d] and the two last edges for |3c[ Note that for each edge of [3c} the 
removal of which doesn't affect the convexity, there is another edge that can be removed 
simultaneously without the loss of convexity, so we have to deal only with the subgraphs 
with 13 edges, whose faces are convex (as otherwise one edge should be added and after that 
another can be removed [or added up to the optimal graph, if the graph is still "non-convex"], 
and we reduce the problem to the one with convex faces by removing an additional edge). 

So, now we have to distinguish a "convex" subgraph from the optimal graph. We use the 
similar technique as for A^ = 6 involving the differential, but we slightly adjust the method 
to have better precision. For every possible subgraph embedding G we find the nearest (in 



the II ■ ||oo-norm on the vertices) graph F isometric to the graphs on Figures 3b -3d (of course 
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we remove the floating vertex from 3b and work on 6 points) and estimate tfie distance. 



Then we use the same estimate on the size of the uniqueness domain for that T (we know its 
coordinates exactly, so we do not need to estimate the differential in the neighborhood). If 
the estimate distance is less than this bound, then we are done, having proved that vertices 
of G coincide with the ones of F. 

Also, using the very original problem statement: to maximize d on the set in ]R^^~^ 
described by inequalities 

' < a;i < 1 
0<yi<l 
dist{pi,pj) > d, 

we've implemented the Karush-Kuhn-Tucker condition checking for proving that there either 
is another edge or the graph is not optimal. It could save us a bit of time, but we had already 
achieved a sufficient precision for the graphs to use the differential estimate by the time these 
check was written. Actually, we used it only for the one graph on 6 vertices. And for the 
another one graph we've had to manually correct the resulting estimates because the basis 
chosen by the program had not been optimal and the better estimates for the coordinates 
could be achieved manually. 

After performing all these steps, all of the graphs have been proved to coincide with 
the optimal ones, except for the only one in the following Figure [6j Here the differential of 
System ([T| degrades and the Karush-Kuhn-Tucker conditions are satisfied, so we have to 
develop an individual strategy for it. 

Let's write down the corresponding system of equations: 

Vlx + V2x + V3x=0 
Vly + V2y + V3y - 1 = 
VAx + V5x + V6^-1 = 

V2x + V^x + Vix + Vjx - nx = 

V2y + V3y + VAy + Vjy -VSy = Q 

V5x + V9x -V7x = 



< V5y + Vgy - Vjy = 

VQx - V3x - 



= 



i "ly 

- V3x - Viix - Viox - ^9. 

- V3y - VUy - VlOy -VQy=0 



(4) 



VGy 

Vl3x - VUx + V2x = 

Vl3y - VUy + V2y = 

Vl2x + Vlx - Vi3^ = 

Vl2y + Vly - Vl3y = 
2 

^ix 



+ v,l-d^ = 



1 < i < 13 



It is a polynomial system of 27 equations in 27 variables. As it was once noticed above, 
there are known methods for isolating the roots of the polynomial systems. However, this 
system also appears to have an infinite number of complex solutions, so these methods fail. 
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Figure 6: The most powerful graph 

if apphed directly. So we try to do some algebraic localization to eliminate the infinite family 
of solutions. Adding the condition (i 7^ 0, rewritten as ci ■ s = 1 by adding an independent 
variable seems to do the thing: now the Maple's Isolate function from RootFinding isolates 
the 203 roots of the modified system. Now, we check that there is a unique root of the system 
even on the set < d < 0.5, and thus it is the one we know, which coincides with the optimal. 
That finishes the proof of the Theorems [2] and |3j 

Remark. One may think that adding such condition d ■ s = 1 could do the thing for all the 
graphs; however, that is far from the general case. 

8 Conclusion 



For iV = 9, we present a result of numerical simulations. 

The code used for solving the problem and its output for = 6 and = 7 is available 
at https://dl.dropbox.coin/u/124992374/index. xhtml. Note that it uses some random- 



ization and has some asserts that may fail sometimes. If the case, the program just needs to 
be restarted for the failed graph. 

Acknowledgment. We wish to thank Alexey Tarasov, Vladislav Volkov and Brittany Fasy 
for some useful comments and remarks, and especially Thom Sulanke for modifying surf tri 
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A Numerical bounds for the graphs 



Here we present the numerical bounds for the edges of Gi and G2 produced by the code 



Go: 



ab: (0.346659 0.34685, -0.200409 -0.200038) 

ac: (0.346668 0.34686, 0.200022 0.200375) 
ad: (-0.265502 -0.265339, 0.299724 0.299871) 
ae: (-0.265475 -0.265366, -0.29985 -0.299742) 
bf: (-0.265466 -0.265357, -0.29988 -0.299747) 
be: (0.387784 0.387866, -0.0997594 -0.0994001) 
cd: (0.38772 0.387911, 0.099223 0.100013) 
cf: (-0.265475 -0.265339, 0.299738 0.299875) 
df: (0.346668 0.34686, 0.200036 0.200367) 



ab: (0.299465 0.30012, 0.265043 0.265794) 

ac: (-0.300166 -0.299428, 0.265001 0.265842) 
ad: (-0.200635 -0.199762, -0.347011 -0.346512) 
be: (0.200062 0.200362, -0.346848 -0.346663) 
be: (0.400342 0.400452, -0.00759171 0.00810797) 
bf: (-0.0997579 -0.0994303, 0.387787 0.387861) 
cd: (0.0992846 0.0998033, 0.387762 0.387948) 
ce: (-0.200335 -0.200062, -0.346848 -0.346687) 



de: (-0.299865 -0.299592, 0.265332 0.265718) 
de: (-0.000109204 0.000109204, 0.4004 0.400427) df: (0.400397 0.400451, -0.00456848 0.00469225) 
ef: (0.346696 0.346832, -0.200325 -0.200083) ef: (-0.299956 -0.299738, -0.265484 -0.265304) 

For Gi we have 5 < 0.002 and e < 0.002. After calculating the inverse of the differential 

we get h > 
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and so the sufficient inequality £ < | holds. 

h 
4- 



Similarly, for G2 we get h> ^ and still £ < 
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B Coordinates of the graphs 



Here are some exact values for the coordinates of the vertices of the optimal graphs: 

_ ^ (1 W 



The following was achieved by writing down the system of equations for the graph as it 
appeared to be (i.e., the edges that seem to be horizontal were assumed horizontal, etc) and 
solving it in MATLAB. 

(0,0), 

(^^^-^^2^3^3 + 2 + ^,-^ (-3x/3 + ^2^373 + 2-1) x/sV 
^73 + ^x727373 + 2,-^ (-3V3 + a/2V3x/3 + 2-i) V3j , 
^,1- ^ V3V2V3V3 + 2 + ^V3 - ^ V2V3V3 + 2^ , 



2 

^ - 4 73 + - v/^75T^, - + - ^ - - V2 x/3 ^ + 2 

1 \ I 5 1^1 



^2 . 4v/3--x/2V3V3 + 2,- + -^--x/2V3^ + 2 

For N = 1 and 8 for all the optimal graphs we just write down the edge coordi- 
nates, as there are only 6 types of them and they look simpler than the point coordinates: 



2i + x/3' 2i + v3rv2i + ^' 21 + ^ r V ' l + \/3y ' Vl + \/3 



22 



